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Abstract 

Generalized convolution symmetries of integrable hierarchies of KP-Toda and 2KP-Toda 
type have the effect of multiplying the Fourier coefficients of the elements of the Hilbert space 
Grassmannian Gr^ + ('H) defined on Ti. = L 2 (S' 1 ) by a specified sequence of constants. The in- 
duced action on the associated fermionic Fock space is diagonal in the standard orthonormal 
base determined by occupation sites and labeled by partitions. The coefficients in the single 
and double Schur function expansions of the r- functions associated to elements W £ Gr^ + (7Y), 
which are its Pliicker coordinates, are multiplied by the corresponding diagonal factors. Apply- 
ing such transformations to matrix integrals, we obtain new matrix models of externally coupled 
type which are also KP-Toda or 2KP-Toda r-functions. More general multiple integral repre- 
sentations of tau functions are similarly obtained, as well as finite determinantal expressions for 
them. 
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1 Introduction: convolution symmetries of r-functions 



Solutions of integrable hierarchies of KP-Toda and 2KP-Toda type are determined by 
their r-functions [23] 121] . KP-Toda r-functions depend on an integer lattice point iV and 
an infinite sequence of abelian flow parameters t = (ti,i2, • • ■), and may be expanded in 
an infinite series of Schur functions s\(t) corresponding to integer partitions A = (Ai > 
A 2 > ••• > 0) 

r(iV,t) = ^7r JV (A)s A (t). (1.1) 

A 

The coefficients 7Tjv(A) are interpretable as Plucker coordinates of the image of an element 
of a Hilbert space Grassmannian [23], [25] consisting of subspaces of the space 7i = L 2 (5' 1 ) 
of square integrable functions on the unit circle in the complex z-plane, embedded via the 
Plucker map in the projectivation of the exterior space AH, (the Fermi Fock space). 

The image of the Grassmannian under the Plucker map consists of decomposable 
elements of ATi, and is determined as the intersection of an infinite number of quadrics cut 
out by the Plucker relations. These are equivalent to an infinite set of bilinear differential 
relations for r(N, t), the Hirota equations [151 123] [21], which are the defining property of 
r-functions. Through the Sato formula for the Baker- Akhiezer function 

<M*, t) = e» u MN ' T M [ *p\ I*- 1 } ■= 2*" 2 , 3." 3 , . . .), (1.2) 

these, extended to include all pairs of points on the integer lattice, are equivalent to the 
KP-Toda hierarchy of differential- difference equations and their associated Lax equations. 

The 2KP-Toda hierarchy [T5], [2Z] can similarly be expressed in terms of r-functions 
depending on N, t and a further infinite sequence of flow parameters t = (ti, t 2 , . . .). 
These admit double Schur function expansions [26] 

r( 2 )(iV, t,t) =Y,52B N (\,n)s x (t)s lt (i), (1.3) 

A fjj 

in which the coefficients B^(X,fi) have a similar interpretation in terms of Plucker co- 
ordinates. They also satisfy an infinite set of bilinear differential Hirota-type relations 
in both sequences of flow variables and difference-differential equations relating different 
lattice points. 

Starting with any given r-function of KP-Toda or 2KP-Toda type, it will be shown in 
the following that new r-functions can be constructed, satisfying the same sets of bilinear 
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relations, having the following Schur function expansions: 



C p (r)(N,t) = J2rx(N)ir N (\)s x (t) (1.4) 

A 

C ( *l(rW)(N,t,i) = ^^r^B^^^t)^), (1.5) 

A /i 

where the factors r\(N), f\(N) are defined in terms of a given pair of infinite sequences 
of no n- vanishing constants {rj} igZ , {rj} i6Z through the formulae 

r x (N) := c r (N) J] r^-i+y, r>(iV) := c f (iV) J] f^,, (1.6) 

(*.j')eA (i,i)£At 

Here the products are over pairs of positive integers E A and E \i that lie 

within the matrix locations represented by the Young diagrams of the partitions A and 
[A, respectively, and 

oo 



Cr(N):=H^, (1.7) 



where 

Pi 



(1.8) 

Pi-i 

The sequence of nonvanishing parameters {p^} may be viewed as Fourier coefficients of a 
function p(z) on the unit circle, or a distribution. It will be shown (Proposition 13. lj) 
that, in terms of the elements of the subspace W C L 2 (S' 1 ) corresponding to a point of 
the Grassmannian, the transformations (II. 4p . (11.51) mean taking a generalized convolution 
product with a function or a distribution p(z) on the unit circle, whose Fourier coefficients 
are related to the ?Ys by: 

n = ^- (1.9) 

Pi-i 

(and similarly for fj). These will therefore be referred to as (generalized) convolution 
symmetries. 

With the usual 2KP-Toda flow parameters (t,t) fixed at some specific values, such 
transformations extend to an infinite abelian group of commuting flows whose parameters 
determine the p^s. This has been used to generate new classes of solutions of integrable 
hierarchies [5], |2H [20] . In the present work, they are studied rather as individual transfor- 
mations, for fixed values of the parameters pi which, when applied to a given KP-Toda 
or 2KP-Toda r-function, produce a new one. Particular cases that implicitly use such 
transformations as symmetries have found applications, e.g., as generating functions for 
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topological invariants related to Riemann surfaces, such as Gromov-Witten invariants and 
Hurwitz numbers [TFJ [18] . 

As an immediate application, we may start with an integral over N x N Hermitian 
matrices: 

Z N (t)= [ dfi(M)e tr ^ tiMi , (1.10) 

where d\i is a suitably defined U (N) conjugation invariant measure on the space H 7Vx7V of 
Hermitian N x N matrices, which is known to be a KP-Toda r-function [29]. Applying the 
convolution symmetry (11.41) for the case where p(z) is taken as the exponential function 
e z on the unit disc, and evaluating at flow parameter values 

t i = ±tr(A i ), t=[A] (1.11) 

for a fixed N x N Hermitian matrix A, we obtain, within a constant multiplicative factor, 
the externally coupled matrix model integral (Proposition 14. ip . 

. N-l 

Z N , ext ([A\) := / dp(M)e trAM = (TTi!) C p (Z N )([A\). (1.12) 

Such integrals arise in a number contexts, such as the Kontsevich-Witten generating 
function [7], the Brezin-Hikami model [61 [301 EI] an d the complex Wishart ensemble 
[H, [28]. More general choices for the function p(z) are shown in Proposition 14.21 to 
also determine KP-Toda r-functions as externally coupled matrix integrals. It is further 
shown, in Proposition 14. 3[ that these matrix model r-functions can be expressed as 
finite N x N determinants. 

Similarly, Hermitian two-matrix integrals with exponential coupling of Itzykson-Zuber 

type ma 

4 2) (t,t)= / d M (M x ) / d[i{M 2 ) e^Z^Mi+uMb+MM) (113) 

are known to be 2KP-Toda r-functions [H dH d21 EH] • Applying the convolution symmetry 
(11. 5p to (11.131) gives an externally coupled two-matrix integral (Proposition [4741 . 

Cfl{zf){[A\, [B\) = [ dv{M x ) [ dp(M 2 ) r r (N, [A], [M 1 ])r,(iV, [B], [M 2 ])e^ M ^\ 

(1.14) 
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where [A] and [B] signify the sequences {itr(A*)} igN + and { jtr(£> i )} ieN + of trace invari- 
ants for the pair of Hermitian matrices A and B and 

r r (N,A,M 1 ) = ^)r A (JV) SA ([A]) SA ([M 1 ]), 

A 

r f (N,B,M 2 ) = J2rx(N)s x ([B])s x ([M 2 ]). (1.15) 

A 

This doubly externally coupled two-matrix model r-function can also be expressed in a 
finite N x N determinantal form (eq. (14.33p . Proposition [4751 . 

This approach can also be extended to more general 2KP-Toda r-functions admitting 
multiple integral representations of the form (14.421) . Applying the convolution symmetry 
(II. 5ft then gives a new 2KP-Toda r-function expressible either as a multiple integral 
(eq. (14.431) . Proposition |4T6|) or as a finite determinant (eq. (14.451) . Proposition |477|) . 

The key to understanding these constructions, and further results following from them, 
is the interpretation of the Sato r-function as a vacuum state expectation value of products 
of exponentials of bilinear combinations of fermionic creation and annihilation operators 
[23| [151 |27] . This well-known construction will be summarized in the next section. 

2 Fermionic construction of r-functions 

We recall here the approach to the construction of r-functions for integrable hierarchies 
of the KP and Toda types due to Sato [231 El], the Kyoto school El [151 EZ] and Segal 
and Wilson [25] . 

2.1 Hilbert space Grassmannian and fermionic Fock space 

We begin with the "first quantized" Hilbert space 7i, which will be identified, as in [25] , 
with the space of square integrable functions on the unit circle 

H = L^S 1 ) =H+ + H- } (2.1) 

decomposed as the direct sum of the subspaces TL+ = span{;z*}j g N and 7Y_ = span{,z _l }j eN + 
consisting of functions that admit holomorphic extensions, respectively, to the interior and 
exterior of the unit circle S* 1 in the complex z-plane, with the latter vanishing at z — oo. 
For consistency with other conventions, the monomial (orthonormal) basis elements of 7i 
will be denoted {e« := z^ 1 ^ 1 }^- 
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Two infinite abelian groups act on 7i by multiplication: 

T + := { 7+ (t) := e^ tiZ % and T_ := { T _(t) := e 2 ^ **"*}, (2.2) 

where t := t 2? • • •) is an infinite sequence of (complex) flow parameters correspond- 
ing to the one-parameter subgroups. More generally, we have the general linear group 
GL(TC) consisting of invertible endomorphisms connected to the identity with well defined 
determinants. (See [25] for more detailed definitions of this and what follows.) 

We consider the Grassmannian Gr-^ + (H) of subspaces W C 7i that are commensurable 
with H + C 7i (in the sense that the orthogonal projection operator n + : W—>TC + is Fred- 
holm while 7r_ : W^Ti- is Hilbert-Schmidt). The connected components of Gr^ + (7i), 
denoted Gt h n(H), N G Z, consist of those W G Gr w+ (7^) for which the Fredholm index 
of 7r + : W^H + is N. These are the GL(TL) orbits of the subspaces 

n% ■= z~ n h+ c n, 

whose elements are denoted W 9; n = g(H+) G Gt h n{TC). The solutions to the KP hier- 
archy are given by the r-function 7jv, fl (t) as defined below, which determines the orbit of 
Wg^ under T + through its Plucker coordinates. The Fermionic Fock space is the exte- 
rior space T := A7i consisting of (a completion of) the span of the semi-infinite wedge 
products: 

\X,N) :=e h Ae i2 A---, (2.3) 

where {ij}je~N+ is a strictly decreasing sequence of integers that saturates, for sufficiently 
large j, to a descending sequence of consecutive integers. This is equivalent to requiring 
that there be an associated pair (A, N) consisting of an integer N and a partition A = 
(Ai, . . . , A<?(a), 0, 0, . . .) of length /(A) and weight |A| = Yli=i where the parts Aj are a 
weakly decreasing sequence of non-negative integers that are positive for i < t(X), and 
zero for i > £(X), such that the sequence is given by 

lj:=Xj -j + N. (2.4) 

In particular, for the trivial partition A = (0), we have the "charge vacuum" vector 

\0,N) =e N -i AeTv-2 A..., (2.5) 

which will henceforth be denoted \N). The full Fock space T thus admits a decomposition 
as an orthogonal direct sum of the subspaces of states with charge A" 

F=®Fn. (2.6) 
yvez 



6 



Denoting by {e*}j e z the basis for H* dual to the monomial basis {e;}i e z for TC, we 
define the Fermi creation and annihilation operators ipi and ip\ on an arbitrary vector 
v £ J 7 by exterior and interior muliplication, respectively: 

ipiV = €i A v, ip\v := ieW, v ETC. (2.7) 

These satisfy the standard canonical anti-commutation relations generating the Clifford 
algebra on TL + TC* with respect to the natural corresponding quadratic form 

[fa ^]+ = Ml = 0, [V-i, = (2-8) 

The basis states |A, N) may be expressed in terms of creation and annihilation operators 
acting upon the charge N vacuum vector as follows [14] 

k 

\X,N) = (-l)Ef=iA H^+^^IN), (2.9) 

i=l 

where (cui, • • -c^klPi, ■ ■ ■ , Ac) is the Frobenius notation (see [16]) for the partition A; i.e., 
ojj is the number of boxes in the corresponding Young diagram to the right of the zth 
diagonal element and A the number below it. 

The Plucker map ^3 : Gr w+ (7i)— >P(jF) takes the subspace W = span(u7 1; w 2 , ■ ■ .) into 
the projectivization of the exterior product of its basis elements: 

*P : span(wi, W2, ■ ■ •) 1— > [w\ A W2 A • ■ ■] , (2-10) 

and may be lifted to a map from the bundle Fr^ + (TC) of frames on Gr^ + (TC) to T 

ty:(w 1 ,w 2 ,...)^w 1 Aw 2 A---. (2.11) 

These interlace the lift of the action of the abelian group T + x TL^TL to Gr w+ (7i) or 
Ft-^ + (TC) with the following representation on T (and its projectivization) 

7+(t) : v i-> 7+(t)f, 7+(t) := e ^ um , v e J 7 , (2.12) 

where 

Hi:=^2MUi, (2-13) 

and t = (ti, t 2 , ■ ■ ■) is the infinite sequence of flow parameters. Similarly, the Plucker maps 
and ^3 interlace the action of the abelian group T_ on Fr^ + (7i) and Gr^ + (7i) with the 
following representation on T 

7_(t) : v i-> 7_(t)f, 7_(t) := e^^ UH ~\ v e J 7 , (2.14) 
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and its projectivization. 

The KP-Toda r-function r g (N,t) corresponding to the element W 9j n G Gr^ + (7i) is 
given, within a nonzero multiplicative constant, by applying the group elements 7+(t) to 
W 9t N, to obtain the T + orbit {W ff) jv(t) := 7+(t)(W ff) jv)}, and taking the linear coordinate 
(within projectivization) of the image under the Plucker map corresponding to projection 
along the basis element \N) 

T g (N,t) = (N\¥(W g , N (t)))- (2-15) 

If the group element g e GL(TC) is interpreted, relative to the monomial basis {ej}j e z, as 
an infinite matrix exponential g = e A of an element of the Lie algebra A G Ql(Tl) with 
matrix elements Aij, then the corresponding representation of GL(7i) on T is given by 

g--^i^ A ^^\\ (2.16) 

where : : denotes normal ordering (i.e. annihilation operators ipj appearing to the right 
when j > and creation operators ipi to the right when i < 0). This gives the following 
expression for Tjv, 9 (t) as a charge iV vacuum state expectation value of a product of 
exponentiated bilinears in the Fermi creation and annihilation operators 

r g (N,t) = (Ntt + (t)g\N). (2.17) 

The equations of the KP and Toda lattice hierarchy are then equivalent to the well-known 
infinite system of Hirota bilinear equations [T5l |23| 121] which, in turn, are just the Plucker 
relations for the decomposable element ^3(W 9i jv(t)). 

Similarly, we may define a 2- Toda sequence of 2-component KP r-functions associated 
to the group element g as follows 

rf (iV, t,t) = (JV|7 + (t)^7_(t)|JV>, (2.18) 
where t = (t\, t 2 , . . .) is a second infinite set of flow parameters. 

Remark 2.1 Note that the image under the Plucker map of the component Gr^N^H.) of the 
Grassmannian Gi-n + (7i) having virtual dimension N is the GL(7i.) orbit of the charged vacuum 
state |JV), and that the GL(TC) action preserves the charge N sectors of the Fock space T . 

2.2 Schur function expansions 

Evaluating the matrix elements of 7+(t) and 7-(t) between the states \N) and |A, N) 
gives the Schur function 

(N\j + {t)\X,N) = (X,N\^{t)\N) = s A (t), (2.19) 
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which is determined through the Jacobi- Trudy formula 

s A (t) = det(/i Ai _, +i (t))| 1 < iij <, (A) (2.20) 
in terms of the complete symmetric functions hi(t), defined by 

oo 

e^ Uzi = ^2hi{t)^. (2.21) 

i=0 

Inserting a sum over a complete set of intermediate states in eqs. (12.171) . (12.181) . we 
obtain the single and double Schur functions expansions 

r fl (JV,t) = ^7r^(A)s A (t), (2.22) 

A 

rf (iV, t,t) = ^ 53s^(A^)a A (t)a M (t). (2.23) 

Here the sum is over all partitions A and /i and 

7r Nt9 (\) = (\,N\g\N) (2.24) 

is the Pliicker coordinate along the basis direction |A, N) in the charge N Fock space Tn 
of the image of the element g(TC+) G Gt h n(?{) under the Pliicker map Similarly, 

B Ni9 (\,n) = (\,N\g\n,N) (2.25) 

may be viewed as the |A,iV) Pliicker coordinate of the image of the element giw^x) G 
Gt h n(TC), where 

w^n := span{e Mi _ i+ Af} G Gt h n(H). (2.26) 

In particular, choosing g to be the identity element I, and using Wick's theorem (or 
equivalently, the Cauchy-Binet identity in semi-infinite form), we obtain [12] 

r?\N,t,i) = (tf|7 + (t)7_(t)|JV) = J> A (t)s A (t) = e^**, (2.27) 

A 

where the last equality is the Cauchy-Littlewood identity (cf. ref. [T6]). 
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3 Convolution symmetries 

3.1 Convolution action on H and Gr^ + (7Y) 

Consider now an infinite sequence of complex numbers {Tj} ig z, and define 

p t := e T % n := — , i G Z. (3.1) 
Pi-i 

In the following, we will assume that the series X^i ^-i converges and that 

lim |rj| = r < 1, (3-2) 



% — >oo 



(although, for some purposes, the latter condition may be weakened). It follows that the 
two series 

oo oo 

P+( z ) = S ^Pi z \ P-( z ) = ^P-i z ~\ (3-3) 

2=0 1=1 

are absolutely convergent in the interior and exterior of the unit circle \z\ = 1, respectively, 
defining analytic functions p±{z) in these regions and that 



R p .= Hp-i (3.4) 



i=l 



converges to a finite value. If the inequality (13.21) is strict, p+{z) extends to the unit 
circle, defining a function in L 2 (S' 1 ). Henceforth, we denote the pair (p + ,p_) by p, where 
the latter can be viewed as a sum p_ + p + in the sense of distributional convolutions, as 
defined below. 

If w G L 2 (5' 1 ) has the Fourier series decomposition 



w{z) = 

l= — oo 



^ ( 3 - 5 ) 



we can define a bounded linear map C p : L 2 (S' 1 )^L 2 (S' 1 ) that has the effect of multiplying 
each Fourier coefficient Wi by the factor p-j-i, and hence each basis element by p». 



oo 



C p (w)(z) = P-i-iWiz\ (3.6) 

i=— oo 

This can be interpreted as applying a pair of generalized convolutions 

C p (w):=(C p+ (w.),C p _(w + )) (3.7) 
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defined by 



C p+ (w^)(z) := lim -V / P+ ((/z)w40d(, (3-8) 

C p _(w + )(z) := lim -L / p-(C/z)w+{0<K, ( 3 - 9 ) 
+o+ 27U2: 7|^| = i +e 



(with the contour integrals taken counterclockwise) to the positive and negative parts of 
the Fourier series 



oo oo 
8=1 1=0 



If p+(z) extends analytically to S 1 , eq. (I3.8P is an ordinary convolution product on the 
circle (in exponential variables). In the examples detailed below, all but a finite number of 
the T_j values vanish for % > 0, and hence the infinite product (13.41) is really finite, p~{z) 
is rational with a pole at z = 1 and the convolution product (I3.9P may be understood on 
S 1 only in the sense of distributions. 

If another pair of functions p± is defined as in (I3.3p . with Fourier coefficients pi then, 
defining the convolution products 

p * p := p_ * p_ + p + * p + 

oo 

p_ * p_( z ) ■= ^2p_ iP _ iZ - 1 

i=l 

oo 

P+ * P+( z ) ■= ^PiPiZ 1 , (3-11) 



i=0 



we have 

Cp*p = CpoC p . (3.12) 



Remark 3.1 Note that the class of generalized convolution mappings defined by (|3.6p - (|3.9|) 
only forms a semi-group since, although they may be invertible, their inverse does not generally 
belong to the same class. It may be extended to a group by dropping the condition (|3.2|) . but this 
will not be needed in the sequel. The linear maps C p : Ti^TC may nevertheless be interpreted 
as elements of GL(TC), and are simply represented in the monomial basis {e{\ by the diagonal 
matrix diagjpj}. They thus belong to the abelian subgroup of GL{7i) consisting of invertible 
elements that are diagonal in the monomial basis. 
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Remark 3.2 Since the Baker- Akhiezer function (jl.2p . evaluated at all values of the param- 
eters t = (tijtg, •••)> spans the shifted element z N {W g ^) £ Gr^o (7i) in the zero virtual 
dimension component of the Grassmannian, the convolution action (13, 8j) . (|3.9h . lifted to the 
Grassmannian, may be obtained by applying its conjugate z N o C p o z _Ar under the shift map 
z _Ar : Gr^o (W)^Gr W Jv('H) to ^Ar(,z,t). But note that, at fixed values of the flow parameters t, 
this does not equal the value of the Baker- Akhiezer function corresponding to the transformed t- 
function as defined below; only the subspaces of 7i that they span, varying over the t values, will 
coincide. This fact will not be used explicitly in the following, but it underlies the geometrical 
meaning of generalized convolutions as symmetries of KP-Toda and 2KP-Toda hierarchies. 

3.2 Convolution action on Fock space 

We now consider the action C x T^*T of the abelian subgroup of GL(7i) consisting 
of diagonal elements in the monomial basis, and associate an element C p G C to each 
sequence {pi}i S z denned as above, such that the Plucker map intertwines the C p action 
with that of C p , lifted to the bundle Frn + (H) of frames over Gr-^ + (7Y), and is equivariant 
with respect to the convolution product ( 13. lip . ( 13.121) and group multiplication in C. 
To do this, we first introduce the abelian algebra generated by the operators 

{ ^ I ^ n (3 ' 13) 
[K^Kj] =0, i,jeZ. (3.14) 

For {pi = e Ti }i<zz as above, define the operator 

d p :=e££-°° TiKi . (3.15) 



Definition 3.1 For each pair (A, N), where N e Z, and X is a partition which, expressed 
in Frobenius notation, is (cvi • • • ctk\Pi • • • Pk), let 

r x (N) := c r (N) J] r N „ i+j = c r (N) (J] -^±2i_ ] , ( 3 .16) 

where 



Pi if iV>0 

c r (N) := { 1 if N = (3.17) 

n=hr if N<0 - 
ii i=N Pi 
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Here the inclusion G A is understood to mean that the matrix location (i,j) corre- 
sponds to a box within the Young diagram of the partition X; i.e., 1 < i < £(X), 1 < j ' < Aj. 
The second equality in A3. 16}) follows from the definition 

It follows that C p acts diagonally in the basis {|A, N)}, with eigenvalues r x (N). 
Lemma 3.1 

C p \X,N)=r x (N)\X,N). (3.18) 



Proof: Since the Fock space basis element |A, N) is an infinite wedge product 

k 

|A, N) = e h A e i2 A • • • = (-l)^*n^Wlr- A -iM. ( 3 - 19 ) 



i=l 



I 



J 



Xj-j + N, (3.20) 



it follows from the definition (12.71) and the normal ordering in (13.131) that the effect of 
the action of e TiKl on |A,iV) is to introduce a multiplicative factor p^ if i > and ej is 
present in the wedge product (13 . 1 9[) or p^ 1 if i < and it is absent, and otherwise no 
factor. Therefore 



A- 



C P \X,N) = C p {-l)^^\[^ N+a ^ N _^_ l \N) (3.21) 



= l^^(f[J^\ \ X ,N 

\U=lP-i 

= c r (N) | IT pN+a * )\X,N) (3.22) 

\l\pN-fk-l) 

= r x (N)\X,N). (3.23) 

Q.E.D. 

Now let W = span{u>i(,2) G L 2 (S' 1 )} ieN + G Gr n+ (H) and view {u'j}j gN + as a frame 

for 

Lemma 3.2 The Pliicker map intertwines the convolution action AS. 6}) and the fr- 
action on JF 

$({C P K) W) = R P C p (${wi} ie -H + ), (3.24) 
with multiplicative factor R p := YliLi P-%- 
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Proof: Applying C p to each element Wi G L 2 (5' 1 ) denning the frame for W G Gr^ + (H) 
just multiplies its Fourier coefficients by the factors p-i-i- It follows that the basis element 
| A, JV) is multiplied by the product of the factors pj. corresponding to the terms it 
contains, as in (I3.18p . Eq. (13.241) then follows from the definition of the Plucker map 
and linearity. Q.E.D. 

Example 3.1 Choose 

oo ,■ 



p + (z) = e 2 = ^-, |z|<l (3.25) 

i=0 

1 00 

p_(z) = -— = Y J Z- i \z\>\, (3.26) 



so 



i' if i < -1, 

I if * ^ 1 
1 if t < 0, 



1 

(nr=i^)(iv); 



where 



is the extended Pochhammer symbol. 
Example 3.2 Choose 



i=0 

and p-(z) again as in (|3.26p . so 



1 if i < 0, 



(3.27) 
(3.28) 



r A (iV)= — y=i „ w ^ if ^(A)<iV (3.29) 



i=i i=i 



(T^ = f>)<^n ICI < !> N<1 (3.31) 



p . = j (a)<S if i > (3>32) 
1 if 2 < —1, 



(3.33) 



rA(JV) = ^ C^'^-^ it <W < A, ,3.34, 



8=0 
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3.3 Convolutions and Schur function expansions of r-functions 

We now consider the KP-Toda tau function 

T Cpg (N,t) = (N\j + (t)C p g\N) (3-35) 

obtained by replacing the group element g in (12.171) by C p g. Such a r-function, obtained 
from r g by applying a convolution symmetry will be denoted 

rc p9 =: C p (r g ). (3.36) 

Introducing a second pair (p + ,p_), defined as in ( 13. 31) . with the Fourier coefficients pi 
replaced by pi, we also consider the 2-Toda tau function 

rgj gC ,(N,t,t) = (N\j + (t)C p gC- p j4t)\N) (3.37) 

obtained by replacing the group element g in (I2.18P by C p gC p , and denote this transformed 
2- To da r-function 

)• (3.38) 

Inserting sums over complete sets of intermediate orthonormal basis states in (13.351) and 
(I3.37p . and defining f\{N) as in (I3.16p . with the factors pi replaced by p~j, we obtain the 
following form for the Schur function expansions ( 12.221) . ( 12.231) . 

Proposition 3.1 The effect of the convolution actions Ii3. 36]) . Ii3. 38\) is to multiply the 
coefficients in the Schur function expansions of Tc p9 (N,t) and Tq § C _(N, t, t) by the diag- 
onal factors r\(N) and r p (N). 

r Cp9 (N,t) = J2rx(N)7r Njg (\)s x (t), (3.39) 

A 

T$ gC _(N,t,i) = $^^r A (iV) J B^(A,/i)f / ,(iV) SA (t) SM (t). (3.40) 

A 

The Pliicker coordinates for the modified Grassmannian elements C p g(?i+) and C p gC p {w p j^) 
are thus 

7T N , Cpg (X)=r x (N)7r N JX) (3.41) 
B NiCpgC - p {\^) = r\(N)B N ,g(\, p)f p (N). (3.42) 
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Proof: This follows immediately from the diagonal form (13.181) of the C action in the 
orthonormal basis {|A, iV)} and the definitions ( 12.241) and (12.251) of the Pliicker coordinates 
7Tjv,c p9 (A) and B NiCpgC ~( A, /i). Q.E.D. 
In particular, setting g — Cp — I, in (I3.40p we obtain 

t®(N, t,t) = £V A (JV) SA (t)s A (t) =: r P (JV,t,t) (3.43) 

A 

where T r (N, t,t) is defined by the second equality. Such r-functions have been studied 
as generalizations of hypergeometric functions in [2TJ [20]. (Cf. also [T2l [13], where the 
notation differs slightly due to the presence of the normalization factor c r (N) in the 
definition (jOBjl of r\(N).) 

In the following, the infinite sequence of parameters t = (t±, t2, ■ ■ ■) will often be chosen 
as the trace invariants of some square matrix M. The sequence so formed will be denoted 



[M\ = <J -tr(M i 



, [M]i := -tr(M l ). (3.44) 
ieN+ 1 



If t and t in (13.431) are replaced by [A] and [B] , respectively, where A and B are a pair of 
diagonal matrices 

A = diag(a 1 ,...,a JV ), S = diag(6i, . . . , 6^) (3.45) 
with distinct eigenvalues, and 

n n 

A(A) := ^ (a, - a,), A(B) := £ (6,, - 6,) (3.46) 

denote the Vandermonde determinants in the variables {a^} and we obtain a simple 
N x N determinantal expression for r r (N, [A], [B]) (cf. (151115] ). 

Lemma 3.3 

T r (N,[A),[B})= E r A (iV) SA ([A]) SA ([5]) (3.47) 

f(A)<Af 



det(p+(a i 6 j )| 1 < ij < iV 



(3.48) 



A(A)A(B) ■ 

Remark 3.3 Although various proofs of this result may be found elsewhere (e.g., cf. [13J), 
we provide a detailed version here, based on the Cauchy-Binet identity in semi-infinite form, 
since it involves some useful further relations. An equivalent way is to use the fermionic form of 
Wick's theorem, which is really just the Cauchy-Binet identity expressed in terms of fermionic 
operators and matrix elements. 



16 



Proof of Lemma l3.3t The Cauchy-Binet identity in semi-infinite form may be expressed 
by considering two iV-dimensional framed subspaces span{Fj}i<j<jv and span{Gi}i<i<N 
of the complex Euclidean vector space £ 2 (N) = span{ej}j G N, identified with 7i + C H = 
L 2 (S' 1 ), by choosing the monomials {z*}jgN as orthonormal basis. The complex inner 
product ( , ) is thus defined by integration 

(F,G):=^-f F(z)G(z- 1 )-. (3.49) 

The Cauchy-Binet identity can then be expressed as 

det(F i ,G j )\ 1 < i j< N = ^2 det(F Xi -i +N j)det(G Xi -i + N,j), (3.50) 

£(\)<N 

where 

F t = J2 F 3i e v G i = J2 ( 3 ' 51 ) 
jez jez 

and the sum is over all partitions A of length £(X) < N, completed so that the N x N 
submatrices F\._ i+N j and Gx^i+Nj are defined by setting Aj = for i > £(X). Since all 
expressions in the sum will be polynomials in the parameters (a iy bj) there is no loss of 
generality in assuming that these lie within the unit disc. We define 

Fi(z) := p+(0iz), G^z) := (1 - b lZ )- 1 (3.52) 

and hence 

F ij = p i (a j ), G lj = (b J ) i . (3.53) 

From the character formula 

det(a X] ~ j+N ) det(b Xj ~ j+N ) 
it follows that the determinant factors on the RHS of (13.501) are 

N 



det(JV^) = det(a?- i+N px^+n) = [Upx^+nJ s x ([A])A([A}), 
det{ gXi - i+Ntj ) = det(b^ i+N ) = s x ([B])A(B). (3.55) 
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From the definitions (I3.16P and (I3.17p . it follows that 




rx(N), (3.56) 

so the RHS of the Cauchy-Binet identity (13.501) is just the RHS of eq. (13.471) multiplied 
by A([A))A([B]). On the other hand, from <KWf . the LHS of (13301 is 

1 f p + (a,iz) dz 



zes 1 z 



= det(p(aj6j)), (3.57) 

which is just the the expression ([MED multiplied by A([A])A([B}). Q.E.D. 
Remark 3.4 Note that, for the case of Example 13.11 eq. (|3.48p becomes the key identity (cf. 

mm) 

1 ,r,n n^-( N Tl\^ d < eaibj) ^<N 



(3.58) 



^(A)<7V 

which, together with the character integral [16] 

d X ,N [ dii H (U)ax([AUXrt]) = s x ([A])s x ([X]), (3.59) 
JueU(N) 

(where <i/u#(C7) is the Haar measure on U(N)), implies the Harish-Chandra-Itzykson-Zuber 
(HCIZ) integral [10] 

Remark 3.5 The condition that the eigenvalues {aj} and of A and B be distinct can be 
eliminated simply by taking limits in which some or all of these are made to coincide. In the re- 
sulting determinantal formulae, like (|3.48p . and those appearing in subsequent sections, in which 
a Vandermonde determinant A(A) or A(B) appears in the denominator, the only modification 
is that the terms in the numerator determinants depending on the a^s and b^s are replaced by 
their derivatives with respect to these parameters, taken to the same degree as the degeneracy 
of their values, while the denominator Vandermonde determinants are correspondingly replaced 
by their lower dimensional analogs. This will not be further developed here, but will be con- 
sidered elsewhere, in connection with correlation kernels for externally coupled matrix models. 
All formulae below in which no Vandermonde determinant factors A(A) or A(B) appear in the 
denominator remain valid in the case of degenerate eigenvalues. 
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4 Applications to matrix models 



We now consider N x N matrix Hermitian integrals that are r-functions, and show how 
the application of convolution symmetries leads to new matrix models of the externally 
coupled type. In the following, let dfi(M), be a measure on the space of N x iV Hermitian 
matrices M e jj ArxAr that is invariant under conjugation by unitary matrices, and such 
that the reduced measure, projected to the space of eigenvalues by integration over the 
group U(N), is a product of N identical measures d/xo on R, times the Jacobian factor 
A 2 (X), 

N 

d/i(UXtf) = TT dfi (x a ) A 2 (X). (4.1) 

'Ueu(N) a=1 

where X = diag(xi, . . . ,Xn). 



J 



4.1 Convolution symmetries, externally coupled Hermitian matrix 
models and r-functions as finite determinants 

It is well-known that Hermitian matrix integrals of the form 

Z N (t)= [ dfi(M)e tr ^ UMl (4.2) 



Men 



N xN 



N r 

= 11/ dix (x a )eZ?=^<A 2 (X), (4.3) 

a=l J * 

are KP-Toda r-functions [29]. The Schur function expansion is 

Z N (t)= J2 ^W^aW, (4.4) 

£(\)<N 

where the coefficients iiN,dfi{^) are expressible as determinants in terms of the matrix of 
moments PU CE21 E3| 

= ( _ 1) |iv(iv-i) iV! det(A^ Ai _ i+i+J v-i)|i<i li <iv (4.6) 

Mij := [ dfi (x)x i+j . (4.7) 
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Now consider the externally coupled matrix model integral (cf. refs. [6| 1281 130| I3T]) 

Z N , ext (A) := / dp{M)e^ AM \ (4.8) 

JMen NxN 

where A e H NxN is a fixed N x N Hermitian matrix. This can be obtained by simply 
applying a convolution symmetry transformation of the type given in Example 13.11 to 
the r-function defined by the matrix integral ( 14.31 ). 

Proposition 4.1 Applying the convolution symmetry C p to the r-function Zjv(t), where 
p+(z) and p-(z) are defined as in VS. 25\) . 113.26}) . and choosing the KP flow parameters 
as t = [A] gives, within a multiplicative constant, the externally coupled matrix integral 

N-l 

C p (Z N )([A}) = (H i^Z N:ext {A). (4.9) 



8=1 



Proof: Using the expansion ([13J) 



e trAX _ d\,N 



£(\)<N \ ^ A 

where 

dx, N = s x (I N ) (4.11) 

is the dimension of the irreducible GL(N) tensor representation of symmetry type A, 
together with the character integral (I3.59P and (14. 5p . it follows that 

2V-1 

Z M (A) = E 7J^^N, d ,(X)s x ([A]) = (l[z\)C p (Z N )\ t=[A] . (4.12) 

Q.E.D. 

More generally, given an arbitrary function p+(z), analytic on the interior of S 1 and 
choosing p-(z) as in (I3.26p . we may define a new externally coupled matrix integral 

Z NtP (A):= I dp(M)r r (N,[AM\), (4.13) 

in which e trAM is replaced by 

r r (N, [M]) := T r (N, [I N ], [M]) = £ d x>N r x (N)s x ([M]). (4.14) 

e(\)<N 
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Then by the same calculation as above, it follows that Zn jP (A) is again just the r-function 
obtained by applying the convolution symmetry C p to Zn, evaluated at the parameter 
values t = [A]. 

Proposition 4.2 Applying the convolution symmetry C p to gives 

C p (Z N ){[A]) = Z NiP (A). (4.15) 

In particular, if we take (p + ,p_) as in Example 13.21 above, we obtain (cf. [13j) 

Z Np (A) = ( TT ^ ] C^ (JV - 1} / dp:(M) det(l - (AM)- a - N +\ (4.16) 

Vf=0 / iMeH^xiV 

showing that this also is a KP-Toda r-function evaluated at parameter values t = [A]. 

Returning to the general finite determinantal formula for Z N ^ P (A) is given by 

the following. 



Proposition 4.3 

^±N(N-l) m 

where 

Gij(p,A):= / dp (x)x l ~ 1 p + (djx). (4.18 



Z »M) = \ f a\ 1 det(G u (p,A))| 1 < jj < 7V , (4.17) 



R 



Proof: Applying the character integral identity (I3.59P to (14.131) gives 

Z NtP (A) = [ dp(M) Y] r x (N)s x ([A])s x ([M}) (4.19) 



1 



d f i (X)A(X)det(p + (a i x j ))\ 1 ^ N (4.20) 
det(G l3 (p,A))\ 1<Li<N , (4.21) 



A(A) 

(_l)i^-i)JV! 



with dj(p,A) defined by (|4.18p . Here, the integration over the U(N) group has been 
performed and Lemma 13.31 has been used in eq. ()4.2Up . Eq. (I4.2ip follows from (14.201) 
by applying the Andreief identity [3j in the form 

N 



JJ / dpo(x m ) I det((/)i(xj)det(il)k(xi) 1<lJ<N = AHdet I / (j)i(x)ipj(x 

V m=l J ) 1<M<JV \J 



l<i,j<N 

(4.22) 
Q.E.D. 
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4.2 Externally coupled two-matrix models 

We now turn to the case of two-matrix models. For simplicity, we only consider the 
Itzykson-Zuber exponential coupling [10], although the same double convolution trans- 
formations may be applied to all the couplings considered in ref. [13]. Using the HCIZ 
identity (13.601) to evaluate the integrals over the unitary groups U(N), we obtain 

4 2) (t,t)= / ^(Mx) / djl(M 2 ) e^UuMi+uM^+M.M,) (4 23) 

= I\ k] I\( f I d ~M e^^< +i ^ + ^A A(X)A(Y), 

k=l a=l J 

where Y = diag(y a , . . . , y^)- This is known to be a 2KP-Toda r-function pfl [21 El 121 131 
[22], with double Schur function expansion 

4 2) (t,t) =^^B Nt ^(\, t j.)s x (t)8 ll (t), (4.24) 

A fi 

where the coefficients B^^dij,{\ p) are N x N determinants of submatrices in terms of 
the matrix of bimoments 

B Ntdtl ^(X,ti) = Y[k\Y[( [ dp (x a ) [ dp Q {y a )e x "A A(X)A(Y)s x ([X])s x ([Y]) 
k =i «=i Jr / 

N 

= (M) JjA;!det(iB Ai _ i+ ^ w _ i+ ^)|i< i j< iV (4.25) 
fe=i 

% := / dfi (x a ) [ djl (ya)e Xaya xy. (4.26) 
Jr Jn 

Now, choosing a pair of elements (p,p), with both p_ and p_ as in (I3.26p . we may 
define a family of externally coupled two-matrix models, by 

Z%l~ p (A, B) := [ dp(M 1 ) [ d~p(M 2 ) r r (N, [A], [M x ])r f (N, [B], [M 2 })e^ M ^\ 

(4.27) 

where A, B are a pair of hermitian N x N matrices. This class may be obtained as the 

~ (2) 

2KP-Toda r-function resulting from applying the convolution symmetry C p ~ p to Z N . 

~ (2) 

Proposition 4.4 Applying the convolution symmetry C p $ to Z N and evaluating at the 
parameter values t = [A], t = [B] gives the externally coupled matrix integral ffJFK l 

C$(ZP)([A]> [B])) = Z { Z rp (A,B). (4.28) 
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Proof: Because of the U(N) x U (N) invariance of the measures dp and dp in (14.270 and 
all factors in the integrand, except for the coupling term e tr ( MlM2 ), we may carry out the 
two U(N) integrations, using the HCIZ identity (13.601) . to obtain a reduced integral over 
the diagonal matrices X = diag(x i . . . , Xjv) , Y = diag(?/i, . . . , y^) of eigenvalues of M 1 
and M 2 , 

4 2) p -(A, B) = H k\ n ( [ dp (x a ) ! dfl (y a ) e^A A(X)A(Y) (4.29) 

k=i «=i V^R / 

xr r (N, [A],[X\)t?(N, [B],[Y]) 
= E E r x( N )BN,d^dfi{^, p)f \{N)s\{[A\)s (t {[B\) (4.30) 

^(A)<iV 1{^)<N 

= C®(2$){[A],[B))). (4.31) 

Q.E.D. 

Since the dependence on A and 5 is U(N) x U(N) conjugation invariant we may 
choose, without loss of generality, A and to be diagonal matrices 

A = diag(oi, . . . ,a N ), B = diag(6 1; . . . , b N ), (4.32) 

We then obtain, as in the one-matrix case, a finite determinantal formula for the 2KP- 
Toda r-function Zjj- JA, B). 

Proposition 4.5 



\ N 

J ^ AB) = ^A(A)A(B) d < Gi ^ * A > B ^ \^i<»> 1 im 

where 



Gij(p,p,A,B) := / dp {x) \ dp (y)e xy p + (a j x)p + (b j y). (4.34) 

JR JR 



Proof 



Z%\~(A,B)= [ dp{M x ) [ dp(M 2 )e^ M ^ 

x E rxWsxQAVaxdMA) £ f^(JV)^([S])« M ([M 2 ]) (4.35) 

^(A)<iV i(n)<N 



A(A)A(S) 
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x det(p + (a fe x/))| 1 < fe Z < Jv det(p + (6 m ?/ n ))| 1 < min < 7V (4.36) 

= A( ( ]) A(B)) det ^^ P> A > B )\i<i,j<N- ( 4 - 37 ) 

In (14.361) . we have used the HCIZ identity ( 13.601) . antisymmetry of the determinants in 
the integrand with respect to permutations in the integration variables (x\, . . . , xn) and 
(yi, ... , yiy) and Lemma 13.31 twice, while in (I4.37p . we have used the Andreief identity 
[3] in the form 

IT / dfi(x m ,y m ) J det(0i(xj)det(^fc(yj) 1<ii<N = AHdet ( / d/i(x, y)4>i{x)^j{y) ) 

m=iJ J i<fci<iv \J / l <ij< N 

(4.38) 
Q.E.D. 

As the simplest example of a 2KP-Toda r-function obtained through Propositions 
14.41 and 14.51 consider the case when the measures d/j, (x) and dp, (y) are both Gaussian, 
and p + and p + are both taken as the exponential function. 
Example 4.1 

dp (x) = e-™ 2 dx, d f i (y) = e-°y 2 dy, p+(x) = e x , p+(y) = ^ '. (4.39) 
Evaluating the Gaussian integrals gives 

Ga = . =e , (4.40) 

3 Vl + 4a 2 V ' 

and hence 

Z N JA) = { ' i '- ii fc=^- e^^^ lK+& ' ) det(e^^). (4.41) 
(l + 4cr 2 )f A(A)A(5) 

The factor e 4 "' 2 - 1 1=1 * * is a linear exponential in terms of the 2KP flow variable ti and 
t<i and hence, through the Sato formula (jl.2p . produces just a gauge factor multiplying the 
Baker- Akhiezer function [25]. Therefore (|4.4ip is just a rescaled, gauge transformed version of 
the 2KP r-function of hyper geometric type appearing in the integrand of the Itzykson-Zuber 
coupled two- matrix model |10j . 



4.3 More general 2KP-Toda r- functions as multiple integrals 

We may extend the above results to more general 2KP-Toda r-functions expressed as 
multiple integrals and finite determinants. To begin with, the following multiple integral 

T$>(N,t,i) = J] ^JJ^dp(x a ,y a )e^^< +i ^ ApT)A(Y), (4.42) 
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where T, f are curves in the complex x- and y-planes and dp(x, y) is a measure on T x T, 
is a 2KP-Toda r-function [13] for a large class of measures dpo(x,y). Applying a double 
convolution symmetry C p - pi with p_ and p_ the same as in (13.261) . gives a new 2KP-Toda 
r-function, also having a multiple integral representation. 

Proposition 4.6 

C%(T$>)(N,t,t) = n (H d ^>y*)) A(X)A(F)r r ( 2 )(iV,t, [X])rf (JV,t, [F]). 

(4.43) 

Proof: This is proved similarly to Proposition 14.41 using the Cauchy-Littlewood iden- 
tity (12.271) twice in the form 

N 

JJ e E£ 1 («.-W.) = £ , A (t), A ([X]) £ S ,(t) S ,([y]). (4.44) 

Q.E.D. 

Evaluating at parameter values t = [A] and t = [B], and applying again Lemma 13.31 
gives the r-function of eq. (14.431) in N x N determinantal form. 

Proposition 4.7 

C { p "l(rS)([A], [B]) = det(G <y (p,p,Ag)li< <J < J y, (4-45) 

(p, p, A, B) := J J dp(x, y)p + (a j x)p + (b j y). (4.46) 



where 



Proof: 

x det (p + (a fc x/ ) ) | j^det (p+ (6 m y n ) ) 1 1 < min < 7V 

JV! 

= A(A)A(B) det(Gij(p, ^ Ajg))l1 ^^ (447) 

where again, we have used the Lemma 13.31 twice and the Andreief identity in the form 
0438]). Q.E.D. 
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This therefore provides a new class of 2KP-Toda r-functions expressible in such a 
finite determinant al form, associated to any pair of curves T, T, together with a measure 
dji on their product, and a pair of functions p+(x) and p+(y), such that the integrals in 
(I4.46P are well defined and convergent. 

Acknowledgements. The authors would like to thank D. Wang for helpful discussions 
relating to this work. 
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